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ABSTRACT 

Spectroscopic measurements of AR NOAA 10871, obtained with the Extreme Ultraviolet Normal 
Incidence Spectrograph (EUNIS) sounding rocket instrument on 2006 April 12, reveal velocity oscil- 
lations in the He II 303. 8A emission line formed at T w 5 x 10 4 K. The oscillations appear to 
arise in a bright active-region loop arcade about 25" wide which crosses the EUNIS slit. The period 
of these transition region oscillations is 26 ± 4 s, coupled with a velocity amplitude of ± 10 km/s, 
detected over 4 complete cycles. Similar oscillations are observed in lines formed at temperatures 
up to T « 4 x 10 5 K, but we find no evidence for the coupling of these velocity oscillations with 
corresponding phenomena in the corona. We interpret the detected oscillations as originating from 
an almost purely adiabatic plasma, and infer that they are generated by the resonant transmission of 
MHD waves through the lower active-region atmospheres. Through use of seismological techniques, 
we establish that the observed velocity oscillations display wave properties most characteristic of 
fast-body global sausage modes. 

Subject headings: Sun: activity — Sun: oscillations — Sun: transition region — Sun: UV radiation 



1. INTRODUCTION 

The discovery of solar oscillations in the 1960s 
(Leighton 1960) and the identification of individual wave 
modes enabled accurate information on the properties 
of the solar interior to be established. These discov- 
eries thereafter developed into the prominent research 
area of helioseismology (Pomerantz et al. 1985; Duvall & 
Harvey 1984; Scherrer et al. 1982). The discovery of 
oscillations in coronal structures raised the prospect of 
using these oscillations as a diagnostic tool to estab- 
lish coronal loop parameters. The field of coronal seis- 
mology subsequently emerged (Roberts 2000; Aschwan- 
den et al. 2004; Nakariakov et al. 2004; Ofman 2007). 
The identification of waves with wavelengths shorter than 
the size of typical coronal structures imposes great de- 
mands on current spectrocscopic and imaging instru- 
ments. The contribution of high frequency waves to the 
heating of coronal structures may be significant (Porter 
et al. 1994). However, current space-borne instruments 
lack the necessary spatial, spectral and temporal resolu- 
tion required to probe such small-scale structures. Ex- 
tensive analyses of photospheric oscillations has shown 
evidence for high-frequency, reconnection-driven acous- 
tic wave trains (Jess et al. 2007b). High-cadence obser- 
vations of the upper solar atmosphere are much more 
scarce, primarily due to telemetry restrictions imposed 
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by space-borne experiment s. On e such study was car- 
ried out by Williams et al. (|2001h who detected a 6 s in- 
tensity oscillation associated with an active region coro- 
nal loop. This oscillation was subsequently interpreted 
as a fast-mode, impulsively generated, magneto-acoustic 
wave (Williams et al. 2002). 

The transition region, with a very dynamic nature over 
an inherently small depth, provides a regime to fully ex- 
plore the potential coupling of underlying wave trains 
to those found in the corona (Doschek 2006). Recent 
work by Gomory et al. (2006J) discusses the detection of 
transition region velocity oscillations with periodicities 
in the range 250-450 s and proposes this a s evid ence of 
downward propagating waves. Doyle et al. ()1998D found 
transition region intensity oscillations with period s rang- 
ing from 200-500 s, while Hansteen et al. (|2000f ) report 
both intensity and velocity oscillations in the transition 
region with periods ranging from 100-300 s. Addition- 
ally, De Pontieu et al. (|2005h find, through use of nu- 
merical models, that oscillations generated in the photo- 
sphere can propagate into the corona providing they are 
guided along an inclined magnetic flux tube. Although 
there have been many detections of transition-region os- 
cillations, space-borne instrumentation lacks the tempo- 
ral resolution required to probe waves with frequencies 
above 33 mHz in this part of the solar atmosphere. The 
reason for this may be the lack of instruments available 
which are sensitive enough to obtain good count statis- 
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Fig. 1. — The top panel shows the EUNIS LW slit position, at the start of the "stare" sequence, after coalignment with eit's 304A image. 
The middle panel shows a snapshot of the LW detector, cropped to only include contributions from the He II emission line, as a function 
of position along the slit. The lower panel plots the integrated He II intensity as a function of position along the slit. Both the middle and 
lower panels are overplotted with vertical dashed lines indicating the boundaries of three bright transition-region features, outlined in § \3\ 
Velocity oscillations were detected in Feature 2. 



tics, while maintaining short exposure times necessary 
for high-cadence observations. In order to probe physi- 
cal processes occurring in the transition region on time 
scales less than 1 min, it is imperative to acquire data 
at the highest possible cadence using a sensitive camera 
system. 

Here we report velocity oscillations detected in the 
transition region using an Extreme Ultraviolet Normal 
Incidence Spectrograph (EUNIS) sounding rocket instru- 
ment developed at NASA's Goddard Space Flight Cen- 
ter. In § 2 we provide a brief outline of the observations, 
while in § 3 we discuss the methodologies used during the 
analysis of the data and the search for reliable oscillatory 
signatures. A discussion of our results in the context of 
velocity oscillations is in § 4, and finally, our concluding 
remarks are given in § 5. 

2. OBSERVATIONS 

EUNIS was launched from White Sands Missile Range, 
New Mexico, at 1810 UT on 2006 April 12. It achieved a 
maximum altitude of 313 km and obtained solar spectra 
and images between 1812 and 1818 UT. EUNIS observed 
NOAA AR 10871 (S07E28) and AR 10870 (S08W02), 
the quiet area between them, and a southern-hemisphere 
coronal bright point within the quiet area. For the 
present study, we are concerned with 75 solar exposures 
taken during a "stare" (fixed pointing) mode near the 



start of the flight, with constant 1.024 s durations. Each 
of the data frames consist of imaged EUV spectra along 
a slit 660" long, together with spectroheliograms (called 
lobe images) formed by 200" wide areas beyond the two 
slit ends which are used for spatial registration. For this 
analysis, we will focus on the long-wave [LW] channel, 
where the spectral coverage ranges from 300A to 370A, 
with an associated cadence of w 2 s. Further observa- 
tions in scanning (spectroheliogram) mode are not dis- 
cussed here. Additional details regarding the EUNIS in- 
strument can be found in Brosius et al. (|2007f ). 

3. DATA ANALYSIS 

Processing the EUNIS images was performed by sub- 
tracting an average dark frame obtained in flight, debi- 
asing each detector row with a linear fit to the unillumi- 
nated ends, correcting for non-linear detector response, 
and applying a flat-field image obtained in the labora- 
tory after the flight. Radiometric calibration, carried 
out post-flight at Rutherford Appleton Laboratory, and 
corrections for atmo spheric extinction, are also included. 
Brosius et al. (2008) discuss validation of the EUNIS rel- 
ative and absolute radiometric calibration. 

Comparison of EUNIS He II 303. 8A lobe images with 
a full Sun He II image obtained at 1812 UT with soho's 
eit reveals that the EUNIS LW slit was tilted about 1.5° 
counterclockwise from the solar E-W axis, and that this 
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Fig. 2. — The slit-averaged wavelength (star symbols) for each 
of the 75 one-second stare exposures, plotted as a function of time 
from the start of the "stare" sequence. For each exposure, the 
slit-averaged wavelength was derived by taking the average of the 
fitted centroid wavelength for each spatial pixel along the slit, as 
outlined in § [3] Notice the decrease in wavelength, with exposure 
number, believed to be caused by thermal effects within the EUNIS 
payload. Overplotted on the diagram, using a solid, bold line, is a 
least-squares fit to the data, which is used to remove this systematic 
wavelength shift. 

angle remained stable in time. The EUNIS LW field of 
view was initially centered at (x = -193", y = -83"), but 
its pointing drifted at a steady rate of about 0.04" s _1 in 
both the x- and y-directions, which corresponds to a net 
drift in each dimension of about 6" during the 167 s spent 
in this mode. Due to the drift in the solar x-direction be- 
ing along the slit axis, we were able to remove this spa- 
tial shift through co- alignment techniques. Implement- 
ing Fourier cross-correlation routines as well as squared 
mean absolute deviations, we were able to provide sub- 
pixel co-alignment accuracy. However, it must be noted 
that sub-pixel image shifting was not implemented due 
to the substantial interpolation errors which may accom- 
pany this technique. Drift of the EUNIS LW field of view 
in the solar y-direction was not compensated for since the 
overall drift is greater than the ~ 2" slit-width. 

After co-alignment of the EUNIS LW spectra, we aver- 
aged over three spatial pixels for a net spatial pixel size 
of 2.8". Line fits were then made using Gaussian profiles 
on linear backgrounds, resulting in values of the centroid 
position, integrated intensity, and FWHM, together with 
their corresponding uncertainties, at each spatial loca- 
tion along the slit and for each EUNIS exposure. Figure 
1 shows the pointing of the EUNIS LW slit near the be- 
ginning of "stare" sequence, along with the simultaneous 
He II slit image and corresponding line intensity as a 
function of position along the slit. Active Region NO A A 
10871, which is the subject of this paper, is located to- 
wards the left-hand side of Figure 1. The EUNIS spectral 
resolution for this first flight was w 200 mA FWHM m 
the LW channel; still, the centroids of strong lines, such 
as with He II reported here, can be determined with a 
statistical uncertainty of ±2 km/s or better (Brosius et 
al. 2007). 

Since EUNIS does not provide an absolute wavelength 
scale, velocities reported here are relative to average val- 
ues measured across individual solar features. For ex- 
ample, Figure [2] shows the slit-averaged He II central 
wavelength as a function of time from the start of the EU- 
NIS stare sequence, and reveals a systematic shift toward 
shorter wavelength. The wavelength drift is modest, less 



than 1 pixel over the whole flight, but real. Its character- 
istics are identical in all observed LW lines, and does not 
depend on the solar feature being observed. These facts, 
coupled with the shape of the drift with time, suggest it 
is due to thermal changes within the payload resulting 
in a tilt of the grating within its dispersion plane. To 
remove this effect, a least-squares linear fit is performed, 
after which the resulting trend is removed from the data 
leaving pure deviations relative to the average value mea- 
sured across the solar feature. Portions of the He II slit 
image were analysed based on their relative magnitudes. 
Figure 1 shows the isolation of three distinct and intense 
features corresponding to a coronal bright point and two 
enhancements in AR NOAA 10871. 

Relative He II velocities as a function of time are gen- 
erated for each of the three features shown in Figure 1. 
Since the spatial averaging of three pixels has already 
been performed to improve the signal-to-noise ratio, no 
additional averaging was required. Thus, spatially re- 
solved velocities are determined for each solar feature, 
allowing for more precise diagnostics. Before temporal 
analysis of the dataset can be implemented, timing cor- 
rections must be considered. During the stare mode in- 
vestigated here, all exposure durations were exactly the 
same, namely 1.024 s. However, time-intervals between 
a small number of exposures varied somewhat. To com- 
pensate for this timing irregularity, all line-fit parameters 
were interpolated onto a constant-cadence time series 
with linear interpolation performed between data points. 
During this interpolation procedure approximately 15% 
of all data points required a timing correction. This pro- 
vides the necessary platform for temporal studies to com- 
mence. 

All spatially-resolved He II velocity measurements were 
passed into Fast Fourier Transform (FFT) and wavelet 
analysis routines. While a FFT searches for periodic 
signatures by decomposing the input signal into infinite 
length sinusoidal wavetrains using a basic exponential 
function, wavelet analysis utilizes a time localised oscil- 
latory function continuous in both frequency and time 
(Bloomfield et al. 2004) and is therefore highly suited in 
the search for transient oscillations. The wavelet chosen 
for this study is known as a Morlet wavelet and is the 
modulation of a sinusoid by a Gaussian envelope (Tor- 
rence & Compo 1998). Strict criteria were implemented 
during wavelet analysis to insure that oscillatory signa- 
tures correspond to real periodicities. The first is a test 
against spurious detections of power that may be due to 
Poisson noise, where the input time series is assumed to 
be normally distributed (consistent with photon noise) 
and following a \ 2 distribution with two degrees of free- 
dom. A 99% confidence level is calculated by multiplying 
the power in the background spectrum by the values of \ 2 
corresponding to the 99th percentile of the distribution 
(Torrence & Compo 1998, Mathioudakis et al. 2003). 

The second criterion applied is a comparison of the in- 
put time series with a large number (1500) of randomized 
time-series with an identical distribution of counts. The 
probability, p, of detecting non-periodic power is calcu- 
lated for the peak power at each timestep by comparing 
the value of power found in the input time series with the 
number of times that the power transform of the random- 
ized series produces a peak of equal or greater power. A 
percentage confidence is attributed to the peak power at 
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Fig. 3. — Wavelet diagram showing the detection of a 26 ± 4 s periodic velocity oscillation in He II, over four complete cycles, originating 
from Feature 2 (see Fig. [TJ- The original He II velocity time series is plotted in a) with arrows indicating the peaks of the oscillation. 
The wavelet power transform along with locations where detected power is at, or above, the 99% confidence level are contained within the 
contours in b). Plot c) shows the summation of the wavelet power transform over time (full line) and the Fast Fourier power spectrum 
(crosses) over time, plotted as a function of period. Both methods have detected a well pronounced 26 s oscillation. The cone of influence 
(COI), cross-hatched area in the plot, defines an area in the wavelet diagram where edge effects become important and as such any 
frequencies inside the COI are disregarded. Periods above the horizontal line (dashed) fall within the COI. 



every time step in the wavelet transform by (1 — p) x 100, 
such that a high value of p implies that there is no peri- 
odic signal in the data, while a low value suggests that 
the detected periodicity is real (see Banerjee et al. 2001). 

Our final wavelet criterion is the exclusion of oscilla- 
tions which last, in duration, less than 1.41 cycles. This 
is consistent with the d ecorrelation time defined by Tor- 
rence k Compo (fl998h . One can distinguish between 
a spike in the data and a harmonious periodic compo- 
nent at the equivalent Fourier frequency by comparing 
the width of a peak in the wavelet power spectrum with 
the decorrelation time. From this, the oscillation lifetime 
at the period of each power maximum is defined as the 
interval of time from when the power supersedes 95% 
significance to when it subsequently dips below 95% sig- 
nificance (McAteer et al. 2004). The lifetime was then 
divided by the period to give a lifetime in terms of com- 
plete cycles (Ireland et al. 1999). Any oscillations which 
last for less than this minimum duration were discarded 
as they may have simply been a spike in the time series. 

Wavelet power, for a variety of oscillatory periods 
which lay above the 95% significance level, was saved 
as the output of this analysis. 

4. RESULTS AND DISCUSSION 

In Figure [3l it is clear to see a strong 26 s oscillatory 
signal, coupled with an amplitude of ±10 km/s, lasting in 
excess of four complete cycles. Through examination of 
the FWHM of corresponding Fourier and global spectra, 
in addition to the boundaries of the 99% confidence-level 
contours, we can associate an error of ± 4 s to the de- 



tected periodic signal. This oscillatory frequency is much 
higher than previously detected and originates from the 
most westerly part of the active region (Feature 2). Sim- 
ilar analysis of the intensity variation of Feature 2 re- 
veals no intensity oscillations at this location. Through 
closer examination of the eit 304A image, the location 
of strong oscillatory power originates from a particularly 
bright part of AR NOAA 10871. Due to the EUNIS 
instrument providing spatially-resolved spectra, pixels of 
the EUNIS slit corresponding to this bright feature (eight 
in total) were analysed in full. Each pixel, corresponding 
to a spatial scale of 2.8" (as outlined in § [3j) , provides high 
oscillatory power at a 26 s periodicity. Furthermore, no 
phase difference between neighboring pixels was found in- 
dicating that this oscillatory phenomena is of substantial 
physical size (in excess of 20"). The fact that the signal 
remains coherent during pointing drifts perpendicular to 
the slit-direction indicates a size of at least 6" in that di- 
mension as well. Additional transition region lines were 
also studied, consisting of the Mg VI 349. 2A (T w 4.0 x 
10 5 K), Ne IV 358.7A (T « 1.6 x 10 5 K) and 
Ne V 359.4A (T sa 3.2 x 10 5 K) emission lines. All 
of these emission lines corroborate the detection of the 
26 s velocity oscillation in Feature 2. However, due to 
these lines being much weaker when compared to He II, 
the resulting reduction in signal-to-noise means that the 
confidence levels associated with these detections are re- 
duced. Nevertheless, a strong 26 s oscillatory signal over 
four complete cycles, with confidence levels exceeding 
85% in each instance, is observed in transition region 
emission lines formed at temperatures up to 4 x 10 5 K. 
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Due to the full-disk nature of eit and MDI images, 
these can be accurately co-aligned by forcing the solar 
limbs to agree (Lopez Fuentes et al. 2006). Through mul- 
tiple coalignments an uncertainty between eit and mdi 
images converges to approximately 0.5 of an eit pixel, or 
1". This is an intermediary step which allows subsequent 
co-alignment of trace images to mdi magnetograms. 
Upon inspection of a co-aligned trace 195A image, 
it is clear that the oscillatory phenomena is gener- 
ated in the immediate vicinity of a coronal loop struc- 
ture (Fig. [4|. This, coupled with the strong underly- 
ing photospheric magnetic- held (top panel of Fig. |4|), 
suggest that this could be an observational signature 
of a propagating MHD wave. Highly magnetic struc- 
tures, such as coronal loops, can act as superior wave 
guides (Edwin & Roberts 1983), which allow the effi- 
cient propagation of MHD waves. To search for up- 
ward propagating waves, we perform the same rigor- 
ous analysis on co-spatial coronal emission lines. For 
this, we choose the strong Mg IX 368. lA (T « 9.5 x 
10 5 K), Fe XIV 334.2A (T » 2.0 x 10 6 K), and 
Fe XVI 335. 4/360. 8 A (T w 2.5 x 10 6 K) emission lines. 
Due to structures at one level in the solar atmosphere not 
always overlying the structures to which they are con- 
nected to at other levels in the atmosphere, we included 
additional EUNIS pixels during our search for coronal os- 
cillations. Due to the EUNIS slit being one-dimensional, 
we are unable to investigate any oscillatory behavior at 
locations North and South of the EUNIS pointing. How- 
ever, we analyse pixels lying to the East and West of the 
detected He II oscillations. We incorporate an extra four 
pixels in each direction into our investigation totalling an 
additional 22.4". After searching for simultaneous oscil- 
lations in these lines, no velocity perturbations could be 
found in the corona. Additionally, no periodic signatures 
in coronal emission intensities could be found implying 
that our detected transition-region oscillations are not 
propagating out into the corona. 

To check the validity, and reliability, of our de- 
tected He II velocity oscillation, we examine the re- 
sults of analysing two additional features mentioned in 
§ [3l One feature corresponds to the central portion of 
AR NOAA 10871 (Feature 1), while the second coin- 
cides with the coronal bright point found to the west 
of the active region (Feature 3). Examining these two 
features in detail reveals no signatures of He II veloc- 
ity oscillations. It is interesting to note that Feature 1 
overlies a region of opposite magnetic polarity to that of 
Feature 2 (Fig. [4j . Additionally, the MDI magnetogram 
shows that the magnetic field is stronger for Feature 2 
than for Feature 1. In accordance with the theory of 
Edwin & Roberts (|1983l ). this may imply that the mag- 
netic structuring around Feature 2 may act as a bet- 
ter waveguide than the magnetic field found surrounding 
Feature 1. Additionally, the velocity fluctuations of the 
He II line are consistent with th e rela tive Doppler veloc- 
ities reported by Brosius et al. (20071). 

Such a scenario leaves open the issue of what causes 
these transi tion r egion velocity oscillations. De Pon- 
tieu et al. (|2003f) state there is a general correspon- 
dence between p-modes and upper transition region os- 
cillations, in duration, perio ds, and locations of oscil- 
latory power. Priest (|l982f) explains that the super- 



FlG. 4. — Top panel shows an MDI magnetogram, the middle 
panel shows an EIT 304A transition region image, and the bottom 
panel shows a TRACE 195A coronal image. Each image was taken 
around the time of the EUNIS launch, and therefore acts as a good 
two-dimensional reference to compare with the one-dimensional, 
spectroscopic results of EUNIS. Ovcrplotted in each diagram are 
the vertical dashed lines corresponding to the boundaries of Fea- 
tures 1 and 2 (§[3] and Fig. [TJ, as well as He II intensity contours 
from the EUNIS LW slit. Notice how Features 1 and 2 are located 
very close to regions of opposite magnetic polarity. 

position of p-mode oscillations usually results in wave- 
trains containing 4 or 5 cycles, after being emitted 
from the p hotosp here. This is corroborated by De Pon- 
tieu et al. (20Q3|) who find transition region oscillations 
typically consisting of wavetrains 3-7 cycles long . The 
cycle li fetim es stated by De Pontieu et al. (2003) and 
Priest (|1982f ) match closely our findings, however, these 
pressure driven modes cannot explain the velocity os- 
cillations presented in this work. Traditionally, transi- 
tion region velocity oscillations, surrounding an active 
region, have been interpreted as the signature of reso- 
nant transmission of chromospheric umbral oscillations 
(Thomas et al. 1987). Previously detected transition re- 
gion velocity oscillations have associated periods of 129- 
173 s (Gurman et al. 1 982) a nd directly corroborate the 
work of Thomas et al. (|1987t) . 

Since there are no low-temperature lines in the EUV 
bandpasses of EUNIS, it provides no direct information 
about chromospheric plasmas. Therefore, we are unable 
to provide a quantitative comparison to the work of Gur- 
man et al. (|1982f ). who detect transition region and chro- 
mospheric velocity oscillation interplay in highly mag- 
netic structures. Nevertheless, the fact that we detect 
pure velocity oscillations, with no intensity variation, im- 
plies that we are dealing with an almost purely adiabatic 
plasma (Rendtel et al. 2003). Models of coupled res- 
onators (Zugzda et al. 1983, 1987; Settele et al. 1999), 
which explain the resonant transmission of MHD waves 
through active-region atmospheres, require regions of 
strong wave reflection, e.g. from the steep tempera- 
ture gradient found within the transitio n reg ion. Con- 
sequently, according to Rendtel et al. (2003), velocity 
oscillations may be observed in the EUV lines of active 
regions as a result of resonant transmission. If intensity 
oscillations were found to accompany transition region 
velocity waves, this would hint at a non-adiabatic sys- 
tem consisting of isothermal waves (Staude et al. 1985). 
Such a distinction is important, as it indicates that the 
solar transition region plasma for Feature 2 is adiabatic 
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in nature, and experiences strong, non-stationary be- 
haviour over the duration of the EUNIS stare mode. 
Therefore, through comp arison of our resul ts to those 
of Gurman et al. dl982D . Rendtel et al. (|2003h and 
Thomas et al. (|1987t ). we can infer that our detected 
transition region velocity oscillations are similar in both 
amplitude and location to those found interacting with 
chromospheric plasma (Brynildsen et al. 1999). We find 
the oscillations propagating in a concentrated magnetic 
field region with much higher frequencies than reported 
to date. 

Of particular importance is the ability to place the ob- 
served He II oscillation onto a dispersion diagram (Ed- 
win & Roberts 1983). From this diagram, it is pos- 
sible to infer the exact mode of propagation in rela- 
tion to fast- or slow-body waves. Several kinds of wave 
modes have been directly observed with the TRACE and 
SOHO spacecraft, as well as with ground-based telescopes 
(Aschwanden et al 1999, Wang et al. 2002). It is of 
particular interest that these wave modes, through use 
of MHD seismology techniques, can be deduced from 
simple plasma parameters such as electron density and 
magnetic-field strength (Nakariakov et al. 1999, Nakari- 
akov & Ofman 2001). One particular mode, the sausage 
mode, has been typically used to interpret short periodic- 
ities of less than one minute. The sausage mode is based 
on an assumption whereby the period of this mode is 
determined by the ratio of the cross-sectional radius to 
the Alfven speed inside the waveguide (Meerson et al. 
1978). The global sausage mode (GSM) is a fast-MHD 
mode which can efficiently modulate the plasma density 
and magnetic field strength. It is also one of the principal 
modes that can be excited in active regions providing the 
magnetic waveguide is sufficiently thick and dense. In or- 
der to assess whether the GSM can be supported in the 
transition region of AR NOAA 10871, we require infor- 
mation on the plasma temperature, electron density and 
magnetic field strength. 

For this assessment, we choose a plasma tempera- 
ture of T k, 5 x 10 4 K (the formation tempera- 
ture of He II), an internal- waveguide electron density of 
n l e w 10 11 5 cm~ 3 , an external- waveguide electron den- 
sity of w 10 10 7 cm -3 and a magnetic field strength 
of 25 G. These values are consistent with previously 
measured transition region values (Parnell et al. 2002, 
Perez et al. 1999). The GSM exists if the dimensionless 
wave number, ka, is greater than the cut-off value, 
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where a is the waveguide radius and j Q — 2.4 is the first 
zero of t he Bes sel function Jo (a;), as detailed in Edwin & 
Roberts £[983). C A0 and C Ae 

arc the internal and exter- 
nal Alfven speeds, while Cso and Cto = CsqCao/ '(do + 
C^o) 1//2 are the internal sound and tube speeds, respec- 
tively. These values can be derived directly from our 
input parameters described above, and fo llowing the cri- 
teria outlined in Nakariakov et al. (|2003h . we can deter- 
mine if this mode can explain the oscillations presented 
here. We derive the associated cut-off value, k c a, to be 
1.12. In order to see whether the observed wave train is a 
fast- or slow-body wave, we relate the detected periodic- 
ity, P, to the internal Alfven speed, whereby a fast-body 
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Fig. 5. — Phase speed u)/k is shown for magneto-acoustic waves 
in a cylindrical waveguide (with radius a) as a function of the 
dimensionless longitudinal wavenumber ka under the transition- 
region/coronal conditions Cab > C*AO > C*T0 > Cso- Solid lines 
show the sausage modes, while the dashed lines indicate kink modes 
(Edwin & Roberts 1983). The solid red box indicates the region 
on the dispersion curve in which the observed He II velocity oscil- 
lations lie. 

wave will be present if: 
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Substituting known values into Equation 2, we find 
that the minimum diameter of the waveguide required 
to channel fast-body waves is s=s 1800 km. We can use 
co-spatial coronal observations to provide an estimate 
of the diameter of the transition-region waveguide. By 
measuring the co-spatial width of coronal loops seen in 
the trace 195A images, we establish an average loop 
width of 7 trace pixels, corresponding to a diameter 
of w 2500 km. This measured width is consistent with 
previous trace statistical loop measurements (Aschwan- 
den et al. 2000, Aschwanden & Nightingale 2005) and is 
greater than the minimum diameter required to channel 
fast-body waves. Accompanying our derived k c a value, 
we deem that a fast-MHD wave (i.e., the global sausage 
mode) could explain the origins of the detected He II ve- 
locity oscillation (FigO. In addition to He II, the warmer 
transition region lines for which periodic velocity oscil- 
lations were detected are also consistent with the global 
sausage mode. Increasing the line-formation tempera- 
ture from 5 x 10 4 to 4.0 x 10 5 K leads to an increase in 
the sound and tube speeds, whereas the internal and ex- 
ternal Alfven speeds remain unchanged. As a result, the 
cut-off value is increased, thus maintaining the validity 
of the GSM as the most favorable mode of propagation. 

However, it must be noted that the cutoff value, k c a, 
is extremely sensitive to the assumed values of magnetic 
field strength and electron density, through their rela- 
tionship with the Alfven speed, Ca- Furthermore, the 
magnetic field strength is almost always an inferred mea- 
sure (Parnell et al. 2002) and exerts a stronger variation 
on the Alfven speed due to a linear dependence between 
them. Therefore, caution must be exercised when evalu- 
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ating the mode of propagation. In this instance, the GSM 
appears to be the most favorable mode of propagation, 
although for a definitive answer, values of magnetic field 
strength and electron density must first be determined 
accurately. 

5. CONCLUDING REMARKS 

We have presented direct evidence of high-frequency 
waves propagating in an active region, at transition re- 
gion heights. We have shown the location of 26 s oscilla- 
tory signatures, with significance levels greater than 95% 
due to the rigorous wavelet criteria enforced in § G3 The 
velocity oscillations detected here originate from a large 
(> 20"), bright He II structure crossing the EUNIS slit on 
the westerly edge of AR NOAA 10871. The highest tem- 
perature at which we observe this oscillation is 4 x 10 5 K, 
but we find no evidence for upward propagation of this 
wave train into the corona, and infer, qualitatively, that 
the detected oscillations could be caused by the reso- 
nant transmission interplay between the chromosphere 
and transition region. We deduce that the observed ve- 
locity oscillations exhibit signatures of a fast-MHD wave 

1 Wavelet software is available at 

http://paos.colorado.edu/research/wavelets/ 



- in particular, that of the fast-body gl obal s ausage wave 
mode. The work of De Pontieu et al. (|2003l ) emphasizes 
the coupling of transition region oscillations to those in 
the underlying photosphere. A natural extension of this 
work is to investigate if such a correspondence also exists 
in the high frequency domain (see Jess et al. 2007a for 
high-frequency photospheric analysis). 
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